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ABSTRACT 

We present mid-infrared spectra of the nearby, Ultraluminous Infrared Galaxy 
NGC 6240 taken with the IRS on Spitzer. The spectrum of NGC 6240 is domi- 
nated by strong fine-structure lines, rotational H2 lines, and PAH emission fea- 
tures. The H2 line fluxes suggest molecular gas at a variety of temperatures. 
A simple two-temperature fit to the S(0) through S(7) lines implies a mass of 
~ 6.7 X lO^Mgat T~ 957K and ~ 1.6 x lO^Mgat T~ 164K, or about 15% of the 
total molecular gas mass in this system. Noteably, we have detected the [NeV] 
14.3//m emission line, with a flux of 5 x 10~^^ erg cm~^ s~^, providing the first 
direct detection of the buried AGN in the mid-infrared. Modelling of the total 
SED from near to far- infrared wavelengths requires the presence of a hot dust 
(T~ 700K) component, which we also associate with the buried AGN. The small 
[NeV] /[Nell] and [NeV]/IR fiux ratios, the relative fraction of hot dust emission, 
as well as the large 6.2/xm PAH EQW, are all consistent with an apparent AGN 
contribution of only 3 — 5% to the bolometric luminosity. However, correcting 
the measured [NeV] fiux by the extinction implied by the silicate optical depth 
and our SED fitting, suggests an intrinsic fractional AGN contribution to the 
bolometric luminosity of ~ 20 — 24% in NGC 6240, which lies within the range 
implied by fits to the hard X-ray spectrum. 
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1. Introduction 

Ultraluminous Infrared Galaxies (ULIRGs), i.e. those galaxies with infrared luminosity 
^IR^ 1O^^L0, have the power output of quasars yet emit nearly all of their energy in the mid 
and far-infrared part of the spectrum. Most ULIRGs are found in interacting and merging 
systems (e.g. Armus, Heckman & Miley 1987; Sanders, et al. 1988; Murphy, et al. 1996), 
where the merger has driven gas and dust towards the remnant nucleus, fueling a massive 
starburst, and either creating or fuehng a nascent AGN (Mihos & Hernquist 1996). 

Spectroscopic observations with the ISO satelhte greatly expanded our understanding of 
the mid-infrared spectra of ULIRGs (e.g., Genzel et al. 1998; Lutz et al. 1999; Rigopoulou 
et al. 1999; Spoon et al. 2002; Sturm et al. 2002; Tran et al. 2001). However, the the 
limitations in sensitivity of the ISO spectrometers, left many ULIRGs, even at relatively low 
redshift, beyond the reach of these methods until now. 

In order to adequately sample the local ULIRG population, we are obtaining mid- 
infrared spectra of a large number of ULIRGs (~ 100) having 0.02 < z < 0.93 with the 
Infrared Spectrograph (IRS) on Spitzer. The entire sample will be observed with the low- 
resolution (~ 80) IRS modules, and about half will be observed with the high-resolution 
{R ~ 650) IRS modules as well. These sources are chosen primarily from the IRAS 1-Jy 
(Kim & Sanders 1998), 2-Jy (Strauss et al. 1992), and the FIRST/IRAS radio-far-IR sample 
of Stanford et al. (2000). In this letter, we present the IRS spectra of NGC 6240, a nearby 
{z = 0.0245), double-nucleus, merging galaxy (Fosbury & Wall 1979; Fried & Schulz 1983; 
Wright, Joseph & Miekle 1984), with an 8 — 1000//m luminosity (Lui as defined by Sanders 
et al 1988) of ~ 7 x lO^^L©. The optical nuclear spectrum of NGC 6240 is classified as 
a LINER (Armus, Heckman & Miley 1989), and the extended optical nebula reveals the 
presence of a starburst-driven superwind (Heckman, Armus & Miley 1987; 1990). X-ray 
observations with ASCA (Turner et al. 1997), Beppo-SAX (Vignati et al. 1999), Chandra 
(Komossa et al. 2003. Ptak et al. 2003), XMM-Newton (Netzer et al. 2005) provide clear 
evidence for the presence of one (or two) AGN behind significant columns of absorbing 
material (Nh = 1 - 2 x 10^^ cm'^). 

Strong [OIV] 25.89;um line emission in the ISO SWS spectrum of NGC 6240 led Lutz et 
al. (2003) to suggest that up to 50% of the infrared energy emitted by NGC 6240 could be 
powered by a buried AGN, yet direct evidence for the AGN in the mid-infrared was lacking 
until now. In this paper we present IRS spectra that show the presence of faint [NeV] 14.3/im 
emission from NGC 6240, and we use fine-structure hne ratios and PAH emission feature 
strengths to estimate the energetic importance of the AGN. Because it is a nearby, merging 
system that has a powerful starburst, a buried (pair of) AGN and a superwind, NGC 6240 
provides us with a close-up view of the critical processes that shape all ULIRGs. Our analysis 
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of the IRS spectra of NGC 6240 will serve to highlight the advantages, and the limitations, 
of IRS spectroscopy for studying luminous, infrared galaxies and buried AGN. 

Throughout the paper, we will adopt a flat, A-dominated Universe {Hq — 70 km s~^ 
Mpc-^ Qm = 0.3, Qa = 0.7). The luminosity distances to NGC 6240 is then 99.2 Mpc, 
assuming z — 0.0245, and one arcsec subtends ~ 450 pc in projection. 

2. Observations 

NGC 6240 was observed with the IRS on 04 March 2004. The IRS is fully described in 
Houck et al. (2004). All four IRS modules were used in Staring mode. Two cycles of 14 sec 
each were obtained for each nod position and each order of the Short-Low (SL) and Long- 
Low (LL) modules. Six cycles of 30 sec each, and four cycles of 60sec each, were obtained 
for each nod position of the Short-High (SH) and Long-High (LH) modules, respectively. A 
high-accuracy blue pcak-up was performed on a nearby offset star before placing the slits on 
NGC 6240. Since the separation of the twin nuclei in NGC 6240 is about 1.5 arcsec, they 
arc both contained within all IRS apertures, the narrowest of which (SL) has a slit width of 
3.6 arcsec. 



3. Data Reduction and Analysis 

All spectra were reduced using the Sll IRS pipeline at the Spitzer Science Center. This 
reduction includes ramp fitting, dark sky subtraction, droop correction, linearity correction, 
and wavelength and flux calibration. One-dimensional spectra were extracted from the two- 
dimensional spectra using the SMART data reduction package (Higdon et al. 2004). The 
one-dimensional spectra are then scaled to a spectrum of a Lacertae, extracted in an identical 
manner using SMART. The SL and LL data have been sky subtracted by differencing the two 
nod positions along the slit, before spectral extraction. As a final step, we have normalized 
the SL and LL ID spectra upwards to match the IRAS 25//m FSC data (Moshir et al. 1990). 
The scale factor was 1.16. The low-resolution IRS spectra are displayed in Fig. la. 

Since the SH and LH slits are too small for on-sht background subtraction, we have 
subtracted the expected background fiux through each sht based on the model of Reach et 
al. (2004). The LH data went through an additional cleaning step to remove residual hot 
pixels, which uses the B-mask (supplied with the BCD data) and the warm pixel mask for 
the given IRS campaign. Hot pixels were interpolated over spatially by using the two nod 
positions, before extraction from 2D to ID spectra. For a full description of this method. 
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see Devest et al. (2005). The SH and LH spectra were then scaled to the corresponding IRS 
low-resolution spectra, using a single scale factor for each module. The scale factors were 
1.15 and 1.07 for SH and LH, respectively. The high-resolution IRS spectra are shown in 
Fig. lb,c. 

4. Results 

4.1. Continuum &c Dust Features 

The IRS low-res spectrum of NGC 6240 shows a steeply rising continuum, heavily 
absorbed at 10 & 18/im by amorphous silicates (Fig. la). There is little or no water ice 
or hydrocarbon absorption from 5 — 7.5/um, as is seen in some nearby ULIRGs - e.g. UGC 
5101 (Armus et al. 2004). Broad PAH (Polycychc Aromatic Hydrocarbon) emission features 
at 6.2, 7.7, 11.3, and 12.6//m are visible. In addition, we detect the much fainter 16.4, and 
17.4//m PAH features, along with the broad "plateau" between 16 — 18yum, which are seen in 
some ULIRGs and nearby starburst galaxies (Sturm et al. 2000, Armus et al. 2004, Smith 
et al. 2004, Brandl et al. 2005). The very weak UIB at 14.22//m (see section 4.2) may also 
be produced by poly-cyclic aromatic hydrocarbons. 

To estimate the contribution of dust emission at different temperatures to the luminosity 
of NGC 6240, and to accurately measure the strengths of the PAH emission features against 
the underlying, silicate-absorbed continuum, we have fit the SED with a mult i- component 
model which includes at least two graphite and silicate dust grain components, PAH emission 
features (fit with Drude profiles), a 3500K blackbody stellar component, and unresolved 
Gaussian emission lines (to fit the fine structure and H2 lines). To extend the SED, we have 
added near-infrared, far-infrared, and sub-mm data to the IRS low-resolution spectra. The 
near-infrared photometry is from Scoville et al. (2000), the 60//m and 100//m data are from 
the IRAS Faint Source Catalog (Moshir et al. 1990), the 120//m, IBOiim, 200//m, ABO/i and 
850^m data are from Klaas et al. (2001), and the 350/im data are from Benford (1999). Our 
fit to the NGC 6240 SED is shown in Fig. 2. For a full description of this model fitting, see 
Marshall et al. (2005). 

Three dust components are required to fit the NGC 6240 continuum spectrum. The 

characteristic temperatures of the three components are 27.1 ± 0.3K (cold), 81.4 ± 1.8K 
(warm), and 680 ± 12K (hot). Each component is composed of graphite and silicate grains 
(with optical properties from Draine & Lee 1984), distributed in size according to the Ry = 
3.1 model of Weingartner & Draine (2001). The hot and warm components are modeled 
as spherically symmetric dust shells which are optically thin in the mid-infrared. The dust 
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grains are distributed with uniform density around the central source, out to a distance of ten 
radial scale factors (where the unit radius is taken to be the distance at which the equilibrium 
temperature of the grains is equal to the component temperature). The temperatures of 
all grains are calculated as a function of their size, composition, and distance from the 
illuminating source. Grain sublimation is taken into account, so that the size-distribution 
can vary with radial distance from the central source, depending upon the temperature of 
the grains. The cold component is modeled as a distribution of grains immersed in a single 
radiation field (i.e., they are not radially distributed from a central source). As such, the 
grains in this component have a range of temperatures due to their varying composition 
and sizes, but not due to their spatial distribution. It is presumed that the cold dust is 
responsible for the absorption of photons from the hot and warm components. Nonetheless, 
the emission from the cold grains is treated as optically thin since they are sufficiently cold 
to radiate predominantly at far-infrared wavelengths, where their opacity is extremely small. 
Note that the range of grain temperatures due to the distribution of grain-sizes in the cold 
component produces an SED that is wider than that obtained with an average grain gray- 
body of the same temperature. In this treatment, there is therefore not a unique temperature 
associated with each component, as there is in a gray-body model. We instead define a 
"characteristic" temperature of a component to be the temperature of the most luminous 
grain size at the distance from the source contributing the majority of the luminosity. This 
luminosity dominating distance corresponds to a tuv ~ 0.5, where approximately half of the 
UV-source photons have been absorbed. Each component therefore contains dust above and 
below the characteristic temperature. With this definition, the characteristic temperature 
roughly corresponds to the expected peak in the dust modified Planck function. 

A hot dust component with a characteristic temperature of about TOOK, as we have 
found in NGC 6240, implies the existence of a significant number of grains near sublimation. 
As suggested by Laurent et al. (2000), an excess of hot dust can be taken as evidence for 
grains heated close to sublimation in the presence of an AGN. In NGC 6240, the presence of 
this hot dust together with the high-ionization lines discussed below, are consistent with this 
interpretation. The hottest grains in thermal equilibrium in the T= 680K dust component, 
if illuminated with an AGN-like spectrum (Sanders et al. 1989) having the bolometric 
luminosity of NGC 6240, would be at a distance of about 1 pc from the central source. The 
hot component accounts for approximately 3.5% of the bolometric luminosity in NGC 6240. 

The optical depths to both the hot and cold dust components are essentially zero, while 
the screen optical depth to the warm component is tw — 5.2 ± 0.1 (effectively the sihcate 
optical depth, since this component dominates at 9.7^m), implying an Ay --^ 95 ± 2 mag. 
This is about a factor of three larger than the lower limit which would be derived from a 
smooth continuum fit anchored at 5.3 — 5.6/xm, lA/im and 34//m, since the un-extinguished 
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model continuum at 9.7//m is dominated by silicate emission. The warm component accounts 
for about 26% of the bolometric luminosity before correcting for extinction, and about 80% 
after correcting for this extinction. 

While a multi-component SED fit is never "unique" , we have compared our best fit to 
alternative models which exclude a hot dust component, and which tie the optical depth to 
the hot dust at the value derived for the warm component. The latter is most similar to 
a spherical shell model wherein the hot dust is closest to the (obscured) nucleus and the 
coldest dust is responsible for most of the obscuration along the line of sight. Both of these 
constrained fits produce significantly poorer results, producing reduced chi-squared values 
that are between 50% and a factor of two larger than our preferred, three component fit 
described above. 

Fits to the multiple PAH emission features in NGC 6240 result in 6.2/7.7 and 11.3/7.7 
flux ratios of 0.21 and 0.25, respectively, suggesting a mix of neutral and ionized grains with 
between 200 — 400 Carbon atoms each (Draine & Li 2001). The 6.2/7.7 flux ratio is between 
that found for NGC 253 and the Circinus galaxy, whereas the 11.3/7.7 flux ratio is about a 
factor of two, larger (Sturm et al. 2000; Draine & Li). The 6.2//m PAH EQW is 0.52, about 
20 — 30% smaller than what is found for pure starburst galaxies (Brandl et al. 2005), but 
much larger (a factor of 5 — 10) than is found for AGN. The PAH-to-L/i? lumninosity ratio 
is about 0.023. In all cases, the extinction to the PAH features is assumed to be zero, as in a 
model where the emission arises from an extended region seen against a (partially) obscured 
continuum. Note, our multiple component continuum fits produce significantly larger PAH 
fiux in features (e.g. the 7.7, 8.6, and 11.3jum features) that are affected most by silicate 
absorptions than do other methods which employ spline fits to the observed continuum, 
even those anchored at relatively "clean" parts of the spectrum. However, the fltting we 
have used to derive the PAH fluxes in Table 1 are directly comparable to the Draine & Li 
(2001) models. This wiU be discussed fully in Marshall et al. (2005). 



4.2. Emission Lines 

The IRS short-high and long-high spectra of NGC 6240 are dominated by unresolved 
atomic, flne-structure lines of Ne, O, Si, S, and Fe, covering a large range in ionization 
potential, (Figs. lb,lc and Table 1). The [Fell] lines at 5.34 and 25.99/xm (partially blended 
with [OIV] 25.89/xm) are particularly strong in NGC 6240 as compared to other ULIRGs 
and starburst galaxies, as has been noted by Lutz et al. (2003). We also detect a faint 
feature at ~ 17.9/im which we attribute to [Fell] emission. Grain destruction, possibly in 
the starburst-driven superwind, is the likely source of the increased gas phase Fe abundance 
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in the cooling ISM in NGC 6240 (Draine & Woods 1990). 

The most striking feature of the NGC 6240 spectrum is the strength of the pure ro- 
tational emission lines of H2. As noted by Draine & Woods (1990), based primarily upon 
the strength of the near-infrared H2 lines, NGC 6240 has both the largest H2 line luminos- 
ity, and the largest ratio of H2 -to- infrared luminosity, of any known galaxy. In the IRS 
spectrum, we detect the S(0) through S(7) lines (see Table 1), except for S(6), which is 
blended with the 6.2/im PAH emission feature. The line flux we measure for the S(0) line 
is 8.6(±2.4) X 10^^"^ erg cm^^ s^^, consistent with the upper limit (< 22 x 10^^"^ erg cm^^ 
s~^) in Lutz ct al. (2003). The H2 excitation level diagram (Fig. 3) provides strong evidence 
for multiple temperature components in NGC 6240, a fact noted previously by Lutz et al. 
Fitting the S(3) through S(7) features, we derive a gas temperature of T~ 957K, and a mass 
of 9.5 X IO^Mq for this warm component. Subtracting this fit off the lower excitation lines, 
we can then fit the S(0) and S(l) lines and derive a gas temperature of T~ 164K, and a 
corresponding mass of 1.6 x IO^Mq. In all cases an ortho-to-para ratio of three is assumed. 
Obviously the coldest component dominates the mass of "warm" H2 in NGC 6240. 

A mass of 1.6 x IO^Mq is about 15% of the molecular gas mass derived from single- 
dish miUimeter CO hue measurements by Solomon et al. (1997). However, it is about 50% 
of the cold molecular gas, or about 1/3 the total (warm plus cold) molecular gas mass, 
within the central 1 kpc as measured by Tacconi et al. (1999), which itself about half of the 
dynamical mass within the central 1 Kpc. The cold gas, as well as the warmer gas traced 
by the near-infrared 1 — S(l) line (Tccza ct al. 2000), is actually centered between the two 
nuclei in NGC 6240, the latter presumably thermally excited via slow shocks triggered by 
either the interaction or the outflowing wind (Tecza et al.. Van der Werf 1996, Egami et al. 
1998). Although we have no spatial information on the H2 emission in our IRS spectrum, 
the de-convolved hnewidths are comparable to those measured by Tecza et al. (500 — 600km 
s~^), suggesting that much of the warm molecular gas we measure is associated with the 
centrally-concentrated gas seen by Tacconi et al. and Tecza et al. between the two nuclei. 

The strongest H2 line in our IRS spectrum, the — S(5) transition at 6.9y[/m, has a 
flux about a factor of 4 — 5 times larger that of the near-infrared 1 — S(l) line (Rieke et al. 
1985; Herbst et al. 1990; Fischer, Smith & Glaccum 1990). The mid-infrared H2 emission 
lines listed in Table 1 alone account for about 0.16% of the total energy output of NGC 6240. 
The model proposed by Draine & Woods (1990) for the strong H2 emission in NGC 6240, 
that of thermal emission from gas heated by X-rays, also predicts an 16.33/im Hj]" rotational 
feature with a flux of ~ 2.5 x 10^^'^erg cm~^ s~^. This feature is unique to the X-ray heating 
model. Unfortunately, this key diagnostic falls near the broad, PAH emission feature at 
16.4//m, which is clearly seen in our IRS spectrum. Since our limits for unresolved features 
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in this part of the spectrum are about 3 x 10~^^erg cm~^ s~^, even without the presence of 
the PAH emission, we cannot confirm or rule out the presence of the hne at the level 
predicted by Draine & Woods. 

Most of the emission hues are resolved in our SH and LH spectra, with the bright 
lines having an average linewidth of ~ 680 km s~^, corrected for instrumental broadening, 
comparable to the linewidths reported by Lutz et al. (2003) from their SWS spectra. There 
is no obvious trend of increased linewidth with wavelength. 

Ratios of the mid-infrared atomic fine structure emission lines can be used to characterize 
the source of the UV photons which ionize the gas and heat the dust (Genzel et al. 1998, 
Lutz et al. 1998, Sturm et al. 2002, Verma et al. 2003). Some features, e.g. the [NcV] 
lines at 14.3 and 24.3/um, imply the presence of an AGN, since the ionization potential (97.1 
eV) is too large to be produced by main sequence stars. The same is not true for [OIV], 
since it takes only 55 eV to ionize O"'""'". While the 25.89yLim is strong in the spectra of AGN, 
it is seen in many starburst galaxies as well (e.g., Lutz et al. 1998, Verma et al. 2003). 
In the IRS spectrum of NGC 6240, [OIV] has a line flux of 27.1 x lO'^^W cm^^ and the 
[OIV] 25.9/[NeII] 12.8 line flux ratio is 0.12. The [SIII] 18.7/[SIII] 33.4 line flux ratio (0.75) 
suggests an electron density of Ug ~ 400cm~^ in the ionized nebula, for T= lO^K. 

An expanded view of the rest-frame 13.7 — 15.0//m region in NGC 6240 (order 14 from 
the SH spectrum) is shown in Fig. 4a. Two broad features are evident - a weak feature 
roughly centered at lA.2fim, and a second, stronger feature, roughly centered at 14.35//m. 
We identify the flrst as 14.22/im PAH emission (Moutou et al. 1996), and the second as 
a blend of the [NcV] 14.32^m and [GUI] 14.36/xm fine structure emission lines. The [NeV] 
line is a direct indication of the AGN nature of NGC 6240. The [GUI] line has been seen in 
the mid-infrared spectra of a number of starburst galaxies (Sturm et al 2000, Spoon et al. 
2000, Devost et al. 2005). The broad, [NeV], [GUI] blend is not seen in other ULIRGs in our 
sample with [NeV] emission (see Fig. 4b and Armus et al. 2005), and is not an instrumental 
artifact. In Fig. 4b we have overlayed the same part of the spectrum from Mrk 273, another 
nearby ULIRG with a strong [NeV] emission line (Armus et al. 2005). In Mrk 273 we see 
neither the 14.22/xm PAH nor the 14.36//m [GUI] line. In NGG 6240 we have flt the [NeV] + 
[GUI] complex with two Gaiissians having fixed positions, and widths equal to that measured 
for the isolated [Nelll] 15.55//.m line (Fig. 4c). Only the relative fiux in the lines is allowed 
to vary. The [NeV] 14.32/xm emission line flux is then 5.1 x 10~^^ W cm~^. After removal of 
the [NeV] and [GUI] lines, the broad feature at 14.2/im is fit with an single, un-constrained 
Gaussian. The [NeV] 14.3/[NeII] 12.8, and the [NeV] 14.3/[NeIII] 15.5 hne fiux ratios are 
then ~ 0.03 and ~ 0.08, respectively. We have not detected the [NeV] 24.32/i,m hne, with 
an upper limit on the line flux of 3.9 x lO^^^W cm-^. The [NeV] 14.3/24.3 line flux ratio is 
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> 1.3, comparable to that found for nearby, bright Seyfert 2 galaxies by Sturm et al. (2002). 

5. Discussion 

We have detected the [NeV] emission line in the IRS spectrum of NGC 6240. This is the 
first reported detection of this feature, and it provides a direct measure of the buried AGN. 
Although NGC 6240 does exhibit large-scale shocks from the superwind, it is unlikely that 
these shocks contribute to the [NeV] emission. As pointed out by Voit (1992), shock models 
tend to enhance the low- ionization lines (e.g., [Nell] or [Sill]). While our measured [Nell] 
12.8/[NeIII] 15.55 line flux ratio of 2.8 is larger than found in most pure AGN models (Voit 
1992), it is much lower than is typically found in shocks, where [NeII]/[Nclll] > 10 (Binette, 
Dopita & Tuohy 1985). It is, however, completely consistent with a moderate excitation 
starburst galaxy (Thornley et al. 2000, Verma et al. 2003). We therefore favor a model 
wherein the [NeV] emission in NGC 6240 comes from a coronal-line region in close proximity 
to one or both of the active nuclei, and that most of the observed [Nell] and [Nelll] emission 
is dominated by the surrounding starburst. 

Pure AGN, and AGN-dominatcd ULIRGs, typically have [NeV]/[NeII] and [OIV]/[NeII] 
line flux ratios of 0.8 — 2 and 1 — 5, respectively (Sturm et al. 2002, Armus et al. 2004). Star- 
burst galaxies, on the other hand typicall have stringent upper limits on [NeV] ([NeV] /[Nell] 
< 0.01), yet a large range in [OIV]/[NeII] flux ratios, from 0.01 - 0.2 (Sturm et al. 2002, 
Verma et al. 2003, Devost et al. 2005). Therefore, both the [NeV] /[Nell] and [OIV]/[NeII] 
line flux ratios in NGC 6240 are significantly below the AGN values, while the [OIV]/[NeII] 
ratio is on the high end, but still consistent with that found in some high-excitation star- 
burst galaxies. If wc use a simple linear mixing model (e.g. Sturm et al. 2002) wherein 
the observed values of [NeV]/[NeII] and [OIV]/[NcII] seen in NGC 6240 are due to "excess" 
[Nell] from the starburst, then both line ratios imply an apparent AGN contribution of about 
3 — 5% to the bolometric luminosity. The detection of a hot dust component, which is not 
seen in IRS spectra of starburst galaxies (Brandl et al. 2005) , and which is responsible for 
~ 3 — 4% of the bolometric luminosity, is additional indirect evidence for the buried AGN. 

Fits to the Beppo-SAX, hard X-ray spectrum (Vignati et al. 1999) suggest an HI 
column density toward the active nucleus of 1 — 2 x 10^^ cm~^, or an Ay ~ 500 — 1000 
mag. Correcting for this HI column implies an intrinsic 2-10 KeV, hard X-ray to bolometric 
luminosity of 0.7 — 2 x 10^^ erg s~^, or a hard X-ray to bolometric luminosity ratio of 
~ 0.025 — 0.075, depending upon the unknown fraction of the X-ray flux which is reflected. 
This range is comparable to other AGN which typically have hard X-ray to bolometric flux 
ratios of 0.01 — 0.2 (see Ptak et al. 2003 for a summary), suggesting that NGC 6240 has 
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a rather "normal" hard X-ray to bolometric luminosity ratio for an AGN (Vignati et al. 
1999). The extinction-corrected, hard X-ray data are therefore consistent with the buried 
AGN producing anywhere from 10 — 100% of the luminosity in NGC 6240, with the larger 
values much more likely. Note, the typical (Elvis et al. 1994) QSO luminosity ratio is about 
0.03. 

While the [NeV]/[NeII] line flux ratio is often treated as reddening independent, this is 
not strictly true. Since the [NeV] is coming from gas directly heated by the AGN, while the 
[Nell] can be dominated by the extended starburst, the distribution of the dust in NGC 6240 
(and other ULIRGs) can have a significant effect on the measured [NcV]/[Ncll] emission- line 
flux ratio. This affect would typically be to lower the observed ratio, causing ULIRGs to 
look more starburst-hke. Could the low [NeV] flux in NGC 6240 be the result of the large 
extinction implied by the X-ray spectrum ? Correcting the measured [NeV] emission by an 
Ay = 1000 mag would produce an intrinsic [NeV]/IR flux ratio in NGC 6240 well above 
any other known AGN or ULIRG, so the hard X-ray derived extinction is not applicable. 
From our SED fit, the extinction to the warm dust is Ay ~ 95 mag. While large, this is 
still likely to be a lower limit to some lines of sight to the circum-nuclear region directly 
excited by the AGN. However, if we assume Ay = 95 mag to the coronal line region, and an 
intrinsic [NeV] /[Nell] line flux ratio of 1.0 for this AGN- heated gas, we derive a corrected 
[NeV] /[Nell] ratio of 0.24 (or 0.20 if we also extinction correct the [Nell] which is excited by 
the AGN). This implies the "intrinsic" AGN contribution is hkely to be between 20 — 24%, 
a number which is within the range estimated from the fits to the extinction-corrected, hard 
X-ray data. 

We would hke to thank Bruce Draine, Aaron Evans, David HoUenbach, Lisa Kewley 
and Bill Reach for many helpful discussions. The comments of an anonymous referee also 
helped to improve the content and presentation of this paper. Support for this work was 
provided by NASA through an award issued by JPL/Caltech. 
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Table 1. Emission Features 



Feature ID 


^vest 


Flux 


EQW 


FWHM 




(/im) 


(10-"erg cm-2s-l) 


(/im) 


(km s-i) 


[Fell] 


5.34 


18.9 (6.9) 


0.024 (0.009) 




Ha S(7) 


5.511 


33.7 (12.5) 


0.042 (0.015) 




PAH 


6.2 


399 (44) 


0.52 (0.06) 




Ha S(5) 


6.909 


95.2 (16.6) 


0.095 (0.016) 




[Aril] 


6.98 


57.0 (15.5) 


0.057 (0.015) 




PAH 


7.7 


1876 (212) 


2.60 (0.28) 




Ha S(4) 


8.025 


36.4 (10.8) 


0.058 (0.017) 




PAH 


8.6 


439 (43) 


1.35 (0.13) 




Ha S(3) 


9.665 


70.9 (1.9) 


0.2393 (0.0041) 


570 (60) 


[SIV] 


10.511 


3.8 (0.3) 


0.0099 (0.0004) 


411 (60) 


PAH 


11.3 


475 (31) 


1.61 (0.10) 




Ha S(2) 


12.279 


39.8 (0.5) 


0.0334 (0.0007) 


542 (60) 


PAH 


12.6 


522 (90) 


0.78 (0.14) 




[Nell] 


12.814 


193.1 (3.7) 


0.1198 (0.0018) 


799 (60) 


PAH 


14.2 


5.8 (1.1) 


0.005 (0.001) 




[NeV] 


14.322 


5.1 (0.9) 


0.0087 (0.0008) 




[ClII] 


14.368 


3.6 (1.0) 


0.0072 (0.0009) 




[Nelll] 


15.555 


70.4 (2.4) 


0.0625 (0.0019) 


828 (60) 


PAH 


16.4 


22.4 (0.8) 


0.0198 (0.0008) 




Ha 8(1) 


17.03 


50.4 (1.5) 


0.0423 (0.0017) 


585 (60) 


PAH 


17.4 


7.3 (1.0) 


0.0055 (0.0010) 




[Fell] 


17.934 


4.6 (0.7) 


0.0039 (0.0006) 


657 (62) 


[SIII] 


18.713 


19.9 (1.4) 


0.0178 (0.0022) 


810 (120) 


[NeV] 


24.318 


< 3.9 


< 0.0021 




[OIV] 


25.890 


27.2 (0.7) 


0.0234 (0.0006) 


674 (60) 


[Fell] 


25.988 


23.1 (0.1) 


0.0210 (0.0003) 


624 (60) 


Ha 8(0) 


28.218 


8.6 (1.2) 


0.0032 (0.0004) 


439 (90) 


[SIII] 


33.481 


26.3 (8.4) 


0.0093 (0.0030) 


537 (60) 


[Sill] 


34.815 


268.9 (10.1) 


0.0981 (0.0005) 


823 (60) 



Note. — ine fluxes and rest-frame equivalent widths for lines as measured in the IRS spectra. Uncertainties in the line fits 
(all single Gaussians except where noted) are listed in parentheses. In all cases, except for the [NeV] 14.32 and [GUI] 14.36 
lines, these uncertainties are the larger of cither the Gaussian fit or the difference in the two nod positions. For [NeV] 14.32 and 
[GUI] 14.36 the quoted errors are conservatively estimated at 20%. The FWHM reported in column five is from a Gaussian fit 
to the SH or LH data, after removal of the instrumental profile, in quadrature. In all cases the error given for the FWHM is the 
larger of either the fit, or the known dispersion due to order curvature and undersampling (~ 60km s~^. FWHM are not listed 
for unresolved emission lines measured only in the SL data (those with wavelengths of A < 9.6/im), or for the PAH features. 
All lines were measured using the SMART spectral reduction package (Higdon et al. 2004), except for the PAH features which 
were measured as a product of the multi-component model fits described in the text. 
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Fig. 1. — IRS Short-Low & Long-Low (top panel), Short-High (middle panel), and Long- 
High spectra (bottom panel) of NGC 6240. Prominent emission lines and absorption bands 
(the latter indicated by horizontal bars) are marked. 
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Fig. 2. — Mult i- component fit to the NGC 6240 spectrum. In the first panel, the fit to 
the entire SED from ~ 1 — 1000//m is shown. In the second panel only the IRS low-res 
wavelength range is shown. In both panels the orange solid line represents the total fit, 
and various spectral components (e.g., stellar light, hot dust, PAH emission, cold dust) 
are highlighted individually. In panel one, only the stellar emission (yellow dotted line) 
and the three continuum dust components, hot (blue, single dot-dashed line), warm (green, 
short dashed line), and cold (red, long dashed line) are shown. In panel two, these four 
components, along with the PAH emission (purple, triple dot-dashed line) are also shown. 
The fitting process is described in the text. Fits to the fine-structure and H2 emission lines 
have been omitted for clarity. 
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Fig. 3. — H2 excitation diagram for NGC 6240, showing for the hnes observed with the 
IRS, the upper level population divided by the level degeneracy, as a function of upper level 
energy. The hnes have not been corrected for extinction. Errors are indicated as vertical 
bars for each line. The fits to the hot (T~ 957K) and cold (T~ 164K) components are 
indicated by a dotted and a dashed line, respectively. 



- 18 - 



-3 



C 
CD 

■o 

X 



0.081 ^ ^ 

: NGC 6240 [NeV] complex 

0.06 - [Nelll] shifted and scaled 
0.04- 
0.02- 

0.00 - — 

-0.02 ' 

: NGC 6240 [NeV] complex 
0.06 - Mrk273 [NeV] shifted and scaled 

0.04- 
0.02- 



(a): 




.> 0.00 



(D 



-0.02 



: NGC 6240 [NeV] complex 
0.06 - 14.22 PAH + [NeV] + [Clll] fit 

0.04- 



0.02 

0.00- 
-0.02 



(b): 





13.6 13.8 14.0 14.2 14.4 14.6 14.8 

rest wavelength (microns) 



15.0 



Fig. 4.— Short-High, order 14 spectrum of NGC 6240, highhghting the PAH, [NeV], [CIH] 
complex at 14.2 — 14.4^m. In (a) we show the NGC 6240 hnes with the [NeHI] 15.55/xm hne 
shifted and scaled to the [NeV] 14.3/im peak. In (b) the NGC 6240 features are compared to 
those seen in another local ULIRG with [NeV] emission, Mrk 273 (see Armus, et al. 2005). 
In (c) we show our fit to the NGC 6240 [NeV] complex with the three Gaussian features 
in dashed lines. In all cases the data have been smoothed with a three-pixel boxcar, and a 
linear fit to the continuum has been subtracted off to place the spectra on a horizontal scale. 



